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1 * 0  INTRODUCTION 

The  Naval  Underwater  Systems  Center  (NUSC)  version  of  the  XY2  Pocential 
Flow  Computer  Program  (XYZPF)  described  in  this  report  was  written  under  an 
ongoing  NUSC  IRIED  project  titled  "Torpedo  Launch  Force  Due  to  Bernoulli 
Effect."  The  main  objective  of  the  project  is  to  create  a  hydrodynamic 
coefficient  computer  code  to  calculate  the  unsteady  forces  and  moments  on  a 
body  travelling  with  five  degrees  of  freedom  (D.O.F.)  motion  (three 
translational  motions,  u,  v,  w,  and  two  rotational,  q,  r)  in  close  proximity 
to  a  wall  or  other  structure. 

The  hydrodynamic  coefficient  code  is  presently  being  developed  and  will 
be  described  in  a  future  technical  memorandum.  Input  to  the  code  will 
require  the  velocity  potential  distribution  on  the  surface  of  the  body  for 
all  five  degrees  of  freedom. 

The  XYZPF  was  chosen  as  the  basis  for  calculating  this  velocity 
potential  distribution  since  it  utilizes  a  very  general  singularity  method 
to  conveniently  calculate  the  potential  flow  about  three-dimensional  bodies 
of  arbitrary  shape.  This  report  describes  the  velocity  potential  code. 

In  particular,  the  DTNSRDC  version  of  the  XYZPF  is  a  computer  program, 
written  in  FORTRAN,  which  utilizes  a  surface  source  density  distribution  on 
the  surface  of  the  body  to  compute  the  irrotational,  incompressible 
potential  flow  about  three-dimensional  arbitrary  shaped  bodies.  The  body 
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surface  is  approximated  by  a  set  of  plane  quadrilaterals  with  a  correspond¬ 
ing  set  of  Fredholm  integral  equations  of  the  second  kind  to  solve  for  the 
source  density  on  each  quadrilateral.  The  equations  are  approximated  by 

(1)  assuming  that  the  source  density  is  constant  in  each  quadrilateral*  find 

(2)  by  applying  the  boundary  condition  (normal  component  of  velocity  is  zero 
for  a  uniform  onset  flow)  at  only  one  point  in  the  quadrilateral,  the 
centroid.  Once  the  set  of  equations  is  solved  for  the  source  density 
distribution,  the  flow  velocity  at  on- and  off-body  points  is  calculated. 
Velocities  are  calculated  for  the  three  unit  translational  onflows  and  any 
other  onflow  specified  by  the  user.  This  code  was  installed  on  a  CDC  6700 
computer. 

The  DTNSRDC  version  of  XYZPF  was  supplied  to  NUSC  by  the  Navy  Engineering 
Software  Systems  (NESS)  and  was  used  as  the  basis  for  the  NUSC  version  of 
the  program.  The  NUSC  version  described  in  this  report  has  been  written  to 
fit  on  a  PDP  11/34  minicomputer  and  has  the  increased  capability  of 
calculating  the  perturbation  and  total  velocity  potentials  at  points  on  and 
off  the  body.  In  addition,  programs  were  created  to  perform  calculations 
for  body  motion  with  two  rotational  degrees  of  freedom,  q  and  r. 

Previously,  difficulties  in  using  the  XYZPF  were  the  tediousness  and 
excessive  time  requirements  involved  in  the  task  of  generating  the  surface 
grid  which  is  used  as  input  to  the  program  to  define  the  surface  of  the 
body.  This  task  has  been  facilitated  by  the  use  of  a  finite  element  mesh 
generator.  The  mesh  generator  used  at  NUSC  has  been  SUPERTAB,  a  product  of 
Structural  Dynamics  Research  Corporation.  An  additional  benefit  of  using 


this  interactive  graphics  mesh  generator  is  the  use  of  its  visual  display  of 
the  grid  as  an  error  check. 


A  post-processor  has  been  written  which  converts  the  node  point 
numbering  scheme  of  SUPERTAB  to  that  required  by  XYZPF.  This  is  described 
further  in  the  "Grid  Generator"  section  of  this  report. 


2,0  THEORETICAL  CONSIDERATIONS 


2.1  VELOCITY  POTENTIAL  AND  BOUNDARY  CONDITIONS 


The  total  velocity  potential  (0)  was  calculated  as  the  sum  of  the  pertur¬ 
bation  potential  (0^)  due  to  the  presence  of  the  body  and  the  potential 
(0^)  due  to  the  onflow: 


where 


4> 

*00 


0  +  d> 

P 


“(V  X  +  V  Y  +  V  Z) 


For  translational  D.O.F. 
For  rotational  D.O.F. 


where  X,  Y,  Z  are  Cartesian  coordinates  of  a  point  in  the 

flow  field. 

^oo  »  *  Vco  are  c^e  tIiree  components  of  flow  field 

X  Y  Z 

velocity  at  infinity. 


The  surface  boundary  condition  for  translational  onflow  is  identical  in 


the  DTNSRDC  and  NUSC  versions  of  XYZPF.  The  program  assumes  the  body  is 
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stationary  and  values  are  computed  separately  for  unit  onflow  in  the 
negative  direction  of  the  X,  Y,  and  Z  axes.  The  surface  boundary  condition 
for  this  case  is  that  the  normal  component  of  velocity  on  the  surface  is 
zero,  i.e.: 


C3£\  m~n  •  V<(>  »  0  TRANSLATION  B.C. 

3")> 

where  n  is  the  unit  outward  normal  vector,  and  V  is  the  gradient. 

For  the  case  of  rotational  motion  of  the  body,  the  onset  flow  is  ro 
and  the  body  is  rotating  with  unit  angular  velocity  in  the  positive 
direction  (right-hand  coordinate  system).  The  surface  boundary  cond 
therefore,  requires  the  normal  velocity  of  flow  on  the  surface  to  equal  the 
normal  velocity  of  the  surface  at  that  point,  i.e.: 


( Mi  *  n  •  7$  *  n  • 

w. 


(u>  xt) 


ROTATIONAL  B.C. 


where  u  is  the  angular  velocity  vector  of  the  body  and  R  is  the 
directional  vector  from  the  center  of  gravity  (C.G.)  of  the  body  to  the 
point  on  the  surface. 

It  should  be  noted  that  the  translation  of  a  body  in  the  positive 
direction  of  a  stationary  flow  field  perturbs  the  fluid  in  the  same  manner  as 
a  stationary  body  in  a  uniform  onflow  in  the  negative  direction.  Therefore, 
perturbation  velocity  potentials  calculated  in  XYZPF  for  uniform  onflow  in 
the  negative  direction  are  identical  to  the  total  velocity  potentials  of  a 
body  moving  in  the  positive  direction  of  a  stationary  flow  field. 
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The  perturbation  velocity  potential  )  is  evaluated  by  one  of 
three  methods  depending  upon  the  ratio  (X/T)  of  the  distance  (X)  between  the 
point  under  consideration  and  the  centroid  of  the  quadrilateral  to  the 
maximum  dimension  (T)  of  the  quadrilateral.  This  approach,  utilizing  three 
ranges,  is  consistent  with  the  remainder  of  the  program  where  integral 
equations  for  the  source  density  were  evaluated  using  the  same  criteria.  As 
stated  in  reference  1,  the  criteria  is:  "If  the  ratio  (X/T)  is  greater  than 
4.0,  the  quadrilateral  is  approximated  by  a  monopole  (as  if  it  were 
concentrated  at  one  point).  If  the  ratio  is  greater  than  2.0  and  less  than 
or  equal  to  4.0,  the  quadrilateral  is  approximated  by  a  quadrupole.  If  the 
ratio  is  less  than  or  equal  to  2.0,  the  integrals  are  evaluated  exactly. 

The  approximate  methods  are  used  because  they  require  much  less  time  than 
the  exact  method. " 

The  perturbation  velocity  potential  equations  for  each  range  were 
derived  in  reference  2.  The  resulting  equations  using  the  terminology  of 
XYZPFP  and  rearranged  to  facilitate  computations  are: 

Exact:  <t>p  =  -  (PR12KCLA1)  ♦  (PR23MCLA2)  +  (PR34)(CLA3) 

+  (PR41 ) (CLA4)  +  (Z)(TVZ) 

Quadrupole:  <t>p  *  (A)(W)  -  (CIXX)(WXX)  +  (CIXYXWXY)  -  (CIYY)(WYY) 

Monopole:  4>p  a  (A)(W) 

W  -  1/  X2  +  Y2  +  Z2  ** 

where:  A  *  Area  of  quadrilateral 
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CIXX,  CIXY,  CIYY  -  Moments  of  inertia  about  the 

origin  of  the  quadrilateral 

CLA1,  CLA2,  CLA3,  CLA4  *  Coefficients  calculated  in  XYZPFP 
PR12,  PR23,  PR34,  PR41  ■  Signed  perpendicular  distance  of 

the  point  from  the  extension  of  a 
side  of  the  quadrilateral  (in 
quadrilateral  coordinate  system). 


wxx 


WXY 

WYY 


32W 
3X3  Y 


X,  Y,  Z  *  location  of  point  in  quadrilateral  coordinate 
system 

TVZ  *  Z  component  of  induced  velocity 
To  validate  the  applicable  range  of  each  formula,  comparisons  were  made 
between  velocity  potentials  computed  by  the  exact  and  approximate  formulas 
for  a  square  quadrilateral  having  unit  source  density. 

Results  are  shown  in  figure  1  where  velocity  potential  is  plotted  versus 
nondimens ionalized  normal  distance  from  the  plane  of  the  quadrilateral. 

With  the  exact  formula  as  a  baseline,  it  is  evident  that,  for  a  square 
quadrilateral,  the  monopole  and  quadrupole  approximations  are  highly 
accurate.  The  two  approximate  methods  converge  to  the  exact  at  the  non- 
dimensional  ized  (X/T)  distance  of  approximately  1.4,  considerably  closer 


Effect  of  Quadrupole  and  Monopole  Approximations  on  Velocity  Potential 
Calculations 
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r. 
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than  the  range  in  which  each  approximation  is  used*  Even  with  these 
excellent  results ,  the  range  criteria  previously  stated  was  not  changed* 

The  reasons  were  to  remain  consistent  with  the  ranges  used  for  other 
calculations  and  to  be  conservative  since  quadrilaterals  of  other  shapes  may 
not  give  rise  to  approximate  values  with  the  same  degree  of  convergence  as 
the  square  quadrilateral* 

The  translation  D*0*F*  programs  automatically  calculate  and  output 
values  of  velocities  and  velocity  potentials  for  the  three  unit  flows  in  the 
negative  direction.  In  addition,  the  user  may  specify  the  vector  of  other 
flows  to  be  considered* 

Similarly,  the  rotational  D*0.F.  programs  also  automatically  calculate 
and  output  values  of  velocities  and  velocity  potentials  for  unit  body  motion  i 

in  the  positive  direction  for  both  rotational  degrees  of  freedom*  However,  1 

no  provision  for  calculations  of  additional  flows  has  been  included* 

Additional  flows  were  intentionally  not  included  since  it  is  a  simple  matter 

to  calculate  them  from  the  unit  flows*  ] 

1 

Since  velocity  potentials  of  each  unit  flow  are  linearly  independent, 
the  total  velocity  potential  for  a  five  degree-of-freedom  flow  becomes: 

$  *  u<b  +  v6  +  wd>  +  q£  +  r$ 
u  v  w  n  q  r 

where  u,  v,  w,  q,  r  are  the  actual  velocity  components 

and  <^,  are  the  unit  velocity  potentials  for 

each  D*0*F. 


S 
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2.2  SYMMETRY  CONSIDERATIONS 

The  translational  D.O.F.  XYZPF  has  the  feature  that  only  nonsym- 
metrical  parts  of  the  surface  have  to  be  defined.  This  feature  capitalizes 
on  the  fact  that  for  uniform  onflows,  geometrically  symmetric  quadrilaterals 
about  a  plane  have  source  densities  equal  in  magnitude.  However,  the  sign 
of  a  reflected  quadrilateral  is  a  function  of  the  plane  about  which  it  is 
reflected  and  also  the  direction  of  motion  of  the  body  or  onflow. 

For  rotational  motion  of  the  body,  the  symmetry  feature  is  only 
applicable  when  the  C.G.  and  centroid  of  the  body  coincide.  This  constraint 
is  imposed  because  this  is  the  only  condition  where  source  densities  of 
geometrically  symmetric  quadrilaterals  are  equal  in  magnitude.  This  is  a 
consequence  of  the  boundary  condition  being  a  function  of  the  displacement 
vector  R  between  the  C.G.  and  the  quadrilateral.  If  the  magnitude  of  R 
varies  between  geometrically  symmetric  quadrilaterals,  then  the  magnitude  of 
the  boundary  condition  will  vary,  resulting  in  unequal  source  densities.  In 
this  case,  the  complete  body  will  have  to  be  defined. 

3.0  GRID  GENERATOR 

As  mentioned  in  the  Introduction,  the  tediousness  and  excessive  time 
requirements  of  manually  generating  the  surface  grid  of  the  body  have  been 
overcome  by  utilizing  a  finite  element  mesh  generator.  Figure  2  is  a 
typical  grid  generated  with  the  use  of  SUPERTAB.  Since  the  body  is 
symmetric,  only  one-quarter  of  the  body  has  been  defined  and  is  represented 
by  those  points  connected  with  lines.  The  remaining  points  are  a  reflection 
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of  the  basic  points  and  are  displayed  only  as  an  aid  to  visualize  the  body 
shape. 

Since  SUPERTAB  has  a  consecutive  node  numbering  scheme  with  only  one 
indice  per  node  and  the  XYZPFP  requires  a  two-indice  per  node  numbering 
scheme,  a  post-processor  called  XYZIN.FTN  was  written  to  accommodate  the 
change.  The  program  is  run  interactively  and  requests  information  as  needed. 

Appendix  A  contains  the  program  listing.  Instructions  on  its  use  are 
contained  in  the  comment  statements  of  the  program  and  are  self  explanatory. 

4.0  PROGRAM  STRUCTURE 

The  NUSC  version  of  XYZPFP  has  been  rewritten  to  fit  on  a  PDP  11/34A 
minicomputer  as  opposed  to  a  mainframe  CDC  6700  computer  for  which 
DTNSRDC*  s  version  was  written. 

The  distribution  of  calculations  into  seven  separate  programs  remains 
the  same  in  both  versions.  The  first  three  programs  calculate  quadrilateral 
geometric  data  and  source  densities.  Program  4  calculates  values  of  velocity 
and  velocity  potential  for  on-bodv  points  while  program  5  performs  similar 
calculations  for  off-body  points.  The  sixth  and  seventh  programs  calculate 
off-body  and  on-body  streamlines,  respectively,  for  translational  motion 
only. 

The  seven  programs  for  the  translational  D.O.F.  are  named  PFP1.FTN 
through  PFP7.FTN.  The  five  programs  for  the  rotational  D.O.F.  are  named 
RPFP1.FTN  through  RPFP5.FTN.  Programs  for  the  translational  D.O.F.  are 
independent  of  the  rotational  D.O.F.  programs  and  must  be  run  separately. 
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In  any  case,  the  first  four  programs  must  be  executed  initially  to 
generate  geometric  parameters  and  calculate  source  densities.  These  data 
are  then  saved  in  three  files  which  are  subsequently  used  when  the  remaining 
programs  are  executed.  The  remaining  programs  can  be  run  in  any  combination 
and  may  be  executed  in  the  same  run  as  the  first  four  programs  or  later 
during  a  separate  run. 

The  files  created  during  execution  of  the  first  four  programs  are  named 
UNIT3.DAT,  UNIT4.DAT  and  POTFILE.DAT  by  the  program.  These  files  are 
necessary  to  execute  programs  5  through  7  and  are  the  default  names  for  the 
files  called  by  these  programs.  In  order  to  manage  these  files,  they  can  be 
appropriately  renamed.  As  discussed  in  the  next  section,  these  non-default 
file  names  can  be  specified  as  the  files  to  be  called  by  the  programs. 

Since  the  translational  and  rotational  D.O.F.  programs  are  separate 
programs,  the  contents  of  the  files  will  be  different.  To  avoid  confusion, 
care  should  be  taken  when  renaming  the  files. 

5 , 0  ORGANIZATION  OF  INPUT  DATA 

Input  data  for  the  NUSC  version  of  XYZPF  is  via  a  computer  file.  The 
input  format  is  basically  the  same  in  both  the  NUSC  and  DTNSRDC  versions, 
therefore,  the  input  description  in  reference  1  is  repeated  here  but  with 
the  changes  required  for  the  NUSC  version. 
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The  format  for  a  new  job  starting  from  program  1  is: 

1st  Line:  RESTART  -  Typed  in  columns  1  to  7 
FLAG  1  -  0  or  1  in  column  19 

0  -  specifies  a  new  job  starting  from  Program  1 
1  -  specifies  a  restart  job  with  Programs  1  through 
4  previously  run 
FLAG  2  -  0  or  1  in  column  21 

0  *  specifies  default  file  names  are  used 
1  -  specifies  non-default  file  names  are  used 
When  Flag  2  on  line  1  is  1  (non-default  file  names)  the  three  renamed  file 
names  are  listed  on  the  next  three  lines  in  the  order: 

UNIT3.DAT  renamed 

UNIT4.DAT  renamed 

POTFILE.DAT  renamed. 

The  next  line  is  for  identif ication.  It  contains  information  to  identify 
the  problem  in  columns  2  to  60. 

The  parameter  line  is  next  and  contains  the  following  information  in  the 
order  given.  The  FORTRAN  format  for  this  line  is  (614,  F8.4,  314). 

NQE  is  the  number  of  quadrilateral  elements  to  be  specified  by 

the  point  lines.  NQE  must  be  in  columns  2,  3,  and  4  and  as 
far  to  the  right  as  possible.  The  maximum  value  permitted 
for  NQE  is  650. 

NSE  is  the  number  of  sections  used.  NSE  must  be  in  columns  6,  7, 

and  8  and  as  far  to  the  right  as  possible. 
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MIX  is  the  maximum  number  of  iterations  to  be  performed  for  the  X 

or  Q  flow.  MIX  must  be  in  columns  10,  11,  12  and  must  be  as 
far  to  the  right  as  possible.  MIX  must  be  less  than  200. 

MIY  is  the  maximum  number  of  iterations  for  the  Y  or  R  flow.  MIY 

must  be  in  columns  14,  IS,  16  and  must  be  as  far  to  the  right 
as  possible.  MIY  must  be  less  than  200. 

MIZ  is  the  maximum  number  of  iterations  for  the  Z  flow.  MIZ  must 

be  in  columns  18,  19,  20  and  must  be  as  far  to  the  right  as 
possible.  MIZ  must  be  less  than  200.  MIZ  is  read  in  during 
a  rotational  D.O.F.  run  but  can  have  any  value  since  it  is 
not  used  during  the  run. 

ISM  i3  the  number  of  planes  of  symmetry. 

ISM  *  0  indicates  there  are  no  planes  of  symmetry. 

ISM  ■  1  indicates  that  Y*0.  is  a  plane  of  synanetry. 

ISM  *  2  indicates  that  Y*0  and  Z*0  are  planes  of  symmetry. 

ISM  *  3  indicates  that  Y*0,  Z*0  and  X*0  are  planes  of  sym¬ 

metry. 

ISM  must  be  in  column  24. 

EPS  is  the  convergence  criteria  used  in  testing  the  convergence 

of  the  iterations.  EPS  must  be  in  columns  25  through  32  and 
should  include  the  decimal  point. 

The  next  line  contains  the  X,  Y,  Z  coordinates  for  the  center  of  gravity. 


The  format  is  3F10.4 


Lines  containing  information  for  the  grid  points  follow.  Each  line  has 
the  format  (3F12.9,  414,  F12.9)  and  contains  the  following  information  for 

one  point  on  the  surface.  See  reference  1  for  grid  requirements. 

XI  is  the  X-coordinate  of  the  point  and  must  be  in  columns  1 

through  12. 

YI  is  the  Y-coordinate  of  the  point  and  must  be  in  columns  13 

through  24. 

ZI  is  the  Z-coordinate  of  the  point  and  must  be  in  columns  25 

through  36. 

NI  is  the  N  index  for  the  point.  NI  must  be  less  than  71  and 

must  be  in  columns  39  and  40. 

MI  is  the  M  index  for  the  point.  MI  must  be  less  than  41  and 

must  be  in  columns  43  and  44. 

NS  is  the  section  identification  number.  NS  may  be  any  positive 

integer  from  1  to  9999  and  must  be  in  columns  45  through  48. 
NS  must  be  different  for  different  sections. 

NE  is  used  to  change  the  direction  of  the  normal  vector.  When 

NE  in  the  first  line  of  a  section  is  not  blank  or  zero,  NI 
and  MI  are  interchanged  for  that  section.  On  the  other 
cards,  NE  is  ignored.  Column  52  is  used  for  NE. 

VN  is  the  normal  component  of  the  velocity  at  the  body  surface. 

Most  problems  have  VN  equal  to  zero.  For  some  problems,  such 
as  inlet  pipes,  it  is  desirable  to  have  a  nonzero  VN  on  some 
sections  of  the  input.  The  value  of  VN  from  the  first  card 
in  a  section  is  used  for  the  entire  section.  Columns  53 


through  64  are  used  for  VN. 
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All  point  lines  for  one  section  must  be  together,  but  within  the  section 
they  may  be  in  any  order. 

Additional  flows  are  calculated  only  for  the  translation  D.O.F. 
progress.  Lines  for  additional  flows  to  be  edited  follow  the  last  point 
line  for  the  last  section.  A  maximum  of  18  flows  may  be  edited.  The  format 
for  these  lines  is  (3F12.9).  On  each  line  are  specified  the  three 
components  of  the  free-stream  velocity  for  one  flow  as  shown  below. 

VXI  is  the  X-component  and  must  be  in  columns  1  through  12. 

VYI  is  the  Y-component  and  must  be  in  columns  13  through  24. 

VZI  is  the  Z-component  and  must  be  in  columns  25  through  36. 

This  is  all  the  information  which  is  required  to  run  the  first  four 
programs. 

The  input  file  format  for  programs  5  through  7  is  identical  to  the 
DTNSRDC  version  except  that  an  additional  line  which  specifies  the  program 
to  be  executed  is  required.  This  line  immediately  precedes  the  standard 
input  data  for  that  program  and  consists  of  the  word  SEGMENTS,  SEGMENT6,  or 
SEGMENT?  (starting  in  column  1)  depending  on  whether  the  program  to  be 
executed  next  is  for  program  5,  program  6,  or  program  7,  respectively.  This 
completes  the  input  data  file  for  a  new  job  starting  from  program  1. 

When  the  job  to  be  run  is  a  restart  job  where  programs  1  through  4  have 
already  been  executed,  the  input  data  file  has  the  structure  of  a  new  job 
data  file  except  the  following  lines  are  excluded: 

Identification  line 
Parameter  line 


Line  with  C.G.  data 


Grid  point  lines 
Additional  flow  lines 
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Appendix  B  lists  two  input  files  for  the  sample  case  of  a  sphere  (the 
same  sample  as  presented  in  reference  1)  with  a  radius  equal  to  1,  a  coarse 
grid  of  12  quadrilaterals  and  using  three  planes  of  symmetry  (see  figure 
3).  The  first  file  is  the  case  where  programs  1  through  5  are  run 
successively,  while  the  second  is  for  the  restart  case  where  only  program  5 
is  run. 


6,0  EXECUTING  THE  PROGRAMS 

To  execute  the  programs,  the  task  (  .TSK)  files  for  each  of  the  PFP  and 
RPFP  programs  must  be  accessible.  In  addition,  the  following  files  must 
also  be  accessible:  XYZ.CMD,  XYZ.TSK,  CREATE. TSK,  and  XYZSUB.OLB.  These 
latter  programs  control  execution  of  the  PFP  and  RPFP  programs  and  contain 
required  subroutines. 

The  procedure  for  executing  a  new  job  or  a  restart  job  is  identical. 

The  procedure  is  simple  and  straightforward  as  evidenced  in  the  sample  of 
figure  4.  The  job  is  run  interactively  on  the  terminal  via  a  command  file. 
To  initiate  the  job,  the  operator  types  "@XYZM  into  the  terminal.  The 
computer  responds  by  asking  whether  the  translational  or  rotational  programs 
are  to  be  executed,  and  the  operator  answers  with  a  "TM  or  "R”  specifying 
translational  or  rotational,  respectively.  The  computer  next  asks  for  the 
input  data  filename  and,  lastly,  the  filename  to  which  the  output  will  be 
written.  This  is  all  that  is  required  of  the  operator. 
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T~tiO  YOU  WANT  TO  RUN  TRANSLATIONAL  Off  ROTATIONAL  XV2  (T  Off  ffi  CS3*  @ 

MifMttiittttittimitftxitmtmitmmmmtmf 


TRANS LAT I ONAL  DECREES  OF  FREEDOff 

XY2  POTENTIAL  FLOW  PROCRAN  -  POP  U/34  VERSION 

PROCRAI1  IS  tiff  mo  TO  $5#  QUADS  ANO  IS#  OFF 
•ODY  POINTS  WITH  50  OFF  I0DV  $TffCA«LIN£S  AND 
50  ON-iOOv  STffEAffUffCS.  PLEASE  REPORT  ANY 
pffottcffs  which  occuff  uith  results  Off  DURING 
EXECUTION.  CAUTION  «  -  THIS  PffOGffAN  USES  A 
LARGE  AffOUNT  OF  OISE  SPACE #  IF  YOU  P«E  GOING 
TO  RUN  A  LARGE  JOI  IT  NIGHT  tE  ADUlSAfLE  TO 
SEE  VOUR  SYSTEffS  PEOPLE  AND  HAVE  A  HANK  DISK 
ASSIGNED  TO  THE  JOt. 


aggga 

isEmeSe.^ 


HARE  CONTAINING  THE  INPUT  DATA 


HARE  TO  UMICH  VOU  UANT  THE  OUTPUT  TO  GO 


TTI 
TT1 
TTI 
TTI 

m 
rri 

TTI 
TTI 
TTI 

>•  <EOF> 

> 


STOP  CffO  OF  XV*  CffCATJON  PHASE 

STOP  END  Of  PFPI 

STOP  END  OF  PFPg 

STOP  END  OF  PFP3 

STOP  END  OF  PFP4 

STOP  END  OF  SEGNENT  5 

STOP  SEGNENT  6  NOT  RUN 

STOP  NO  SCRATCH  FILES  CREATED  FOR  SECTION  7 
STOP  SEGNENT  7  NOT  RUN 


NOTE:  ITEM  IN  SLOCKS  ARE  INPUT  BY  OPERATOR 


Figure  4.  Program  Execution  Procedure 
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If  a  program  is  not  being  executed,  the  computer  acknowledges  this  and 
notifies  the  operator  via  the  terminal.  If  a  program  is  executed,  the 
computer  makes  acknowledgement  upon  completion  via  the  terminal. 

A  file  called  FILENAMES.DAT,  created  during  the  execution  of  the 
programs,  remains  in  the  user*s  directory  at  the  end  of  the  job.  This  file 
is  useful  only  during  the  execution  of  the  immediate  job  and  may  be  deleted 
upon  completion  of  the  job.  A  new  file  is  automatically  created  at  the 
beginning  of  each  new  and  restart  job. 

7 . 0  OUTPUT 

The  sample  case  of  the  sphere  with  a  radius  equal  to  1  was  run  for  both 
the  translational  and  rotational  XYZPF  programs.  The  output  listings  are 
presented  in  appendix  C. 

Output  from  the  translational  programs  are  similar  to  the  DTNSRDC 
version  except  for  the  addition  of  total  and  perturbation  velocity 
potentials  for  each  on-  and  off-body  point.  A  similar  output  listing  is 
also  generated  for  the  rotational  D.O.F.  calculations.  However,  since  the 
total  velocity  potential  equals  the  perturbation  velocity  potential,  only 
one  value  is  given.  The  output  table  from  the  rotational  programs  also 
differs  in  that  pressure  coefficient  (CP)  is  not  given  for  the  rotational 
D.O.F.  This  is  because  pressure  coefficient,  as  calculated  in  the 
translational  D.O.F.  programs,  is  not  applicable  to  the  rotational  case. 
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8.0  RESULTS 

Results  of  velocity  potential  calculations  for  the  sample  case  of  the 
sphere  in  a  uniform  stream  is  presented  in  figure  5.  It  is  evident  that 
even  though  this  sample  has  a  coarse  grid,  excellent  results  were  obtained. 

A  second  sample  was  run  using  a  2:1  ellipsoid.  Since  this  body  is 
symmetric  about  three  planes,  only  one-eighth  of  the  body  was  modeled.  A 
coarse  grid  of  only  32  quadrilaterals  was  chosen. 

Figures  6  through  8  present  results  for  two  translational  onflows  and 
one  rotational  motion  of  the  ellipsoid.  Again,  as  with  the  sphere  example, 
excellent  results  were  obtained. 
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APPENDIX  C 

OUTPUT  LISTING  -  SPHERE  SAMPLE 
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